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Abstract 

We investigate freeze-out in hydrodynamic models for relativistic heavy-ion colli- 
sions. In particular, instantaneous freeze-out across a hypersurface of constant tem- 
perature ("isothermal" freeze-out) is compared with that across a hypersurface at 
constant time in the center-of-momentum frame ("isochronous" freeze-out). For one- 
dimensional (longitudinal) expansion the rapidity distributions are shown to differ 
significantly in the two scenarios, while the transverse momentum spectra are remark- 
ably similar. We also investigate the rapidity distribution in greater detail and show 
that the Gaussian-like shape of this distribution commonly associated with the Landau 
expansion model in general emerges only if one neglects contributions from time-like 
parts of the isothermal freeze-out hypersurface. 
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1 Introduction 



Relativistic hydrodynamic models are frequently applied to study high-energy heavy-ion 
collisions |jl|, 0, because they simply represent (local) energy-momentum and charge conser- 
vation. Their main advantage is that they provide a direct link between space-time dynamics, 
i.e., observable flow phenomena, and the nuclear matter equation of state. They are there- 
fore a unique tool to investigate e.g. effects of phase transitions, such as the deconfinement 
transition to the quark-gluon plasma, on the dynamics of the reaction ^ |, ||. 

Relativistic hydrodynamics requires the assumption of local thermodynamical equilibri- 
utttQ i.e., it describes the evolution of an ideal fluid 0. This assumption obviously breaks 
down in the very early stages of ultrarelativistic heavy-ion collisions (the stopping power of 
nuclear matter is insufficient to produce immediate thermalization of matter in the collision 
zone) and in the very late stages of heavy-ion collisions at all energies. In this stage, the 
interaction rate between individual particles becomes too small, or equivalently, their mean 
free path between collisions becomes too long, to maintain local thermodynamical equilibri- 
um in the expanding system, and non-equilibrium processes become dominant. Eventually, 
particles will decouple completely and stream freely until they reach the detector. 

In order to calculate experimental observables it is essential to incorporate this so-called 
freeze-out into the hydrodynamical simulation of a heavy-ion collision. The long mean 
free path or the small interaction rate, respectively, suggest a treatment of the freeze-out 
via kinetic theory [|^. Up to now, however, no attempt has been made to combine the 
hydrodynamical and the kinetic description to model the late expansion stage of a heavy- 
ion collision. The most simple and therefore often used method is to assume an instantaneous 
transition from fluid-dynamical motion to free-streaming. Note that microscopic models for 
heavy-ion collisions |^ do not encounter this problem, since they follow the trajectories of 
individual particles. On the other hand, their disadvantage is that they do not provide a 
direct link to the nuclear matter equation of state and, when based on binary scatterings, 
become invalid for higher densities. 

In this paper we present a systematic comparison of two different instantaneous freeze-out 
scenarios presently used in the literature, namely the freeze-out across a space-time hyper- 
surface defined by a constant temperature (referred to as "isothermal" freeze-out in the 
following) and that across one defined by a constant time in the center-of- momentum (CM) 
frame of the system, respectively (subsequently termed "isochronous" freeze-out). Contri- 
butions from time-like and space-like parts of the hypersurfaceQ are discussed separately and 
it is shown that the often used approximation of neglecting the time-like contributions in 
the isothermal freeze-out scenario leads to serious distortions of the rapidity spectra, and 
moreover violates the conservation laws. 

The above freeze-out scenarios are not strictly self-consistent, because the hypersurface 
is determined after the fluid-dynamical problem has been solved in the whole forward light 
cone. In other words, the influence of the freeze-out of particles from the (time-like or 

"'^A well-defined theory of relativistic dissipative hydrodynamics, which would allow for deviations from 
local thermodynamical equilibrium, does not yet exist, see for instance the discussion in 
definition of time-like and space-like surfaces is given in Section 3. 
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space-like) surface of the fluid on the dynamical evolution is neglected. Although attempts 
to treat this problem analytically have been made 0, up to now no generally accepted 
solution exists. 

This paper is organized as follows. In Section 2 we present the underlying relativis- 
tic hydrodynamic expansion model. In Section 3 we introduce the two different freeze-out 
scenarios. Section 4 contains a discussion of single-inclusive rapidity distributions and trans- 
verse momentum spectra. In Section 5 we present our results and close with a summary in 
Section 6. We use natural units, Ti = ks = c = 1. 



2 The expansion model 

Since we want to focus on the effects of different freeze-out scenarios on the final particle 
spectra, we choose a particularly simple and transparent hydrodynamic expansion model, 
namely the one-dimensional expansion of a finite slab of matter with an ultrarelativistic ideal 



gas equation of state. This model was first studied by Landau et al. []TOl- The expansion 
with Landau initial conditions can to first approximation serve to describe the evolution 
of (net) baryon-free matter in the central region of ultrarelativistic heavy-ion collisions 
(v^ ~ lOOAGeV). In the following we briefiy outline the essential ingredients and results 
of the Landau expansion model. 

In general, relativistic hydrodynamics implies (local) energy-momentum conservation, 

d^T^^" = , (1) 

where for an ideal fiuid the energy-momentum tensor reads 

T'^- = (e + p)u^'u'' - pg>"' . (2) 

Here, e and p denote the energy density and pressure in the local rest frame of matter, 
gfif _ ciiag(+l, —1,-1, —1) is the metric tensor, and u'^ is the fluid 4-velocity {u^u^ = 1). 
The system of eqs. (|l|, |^) is closed by choosing an equation of state for the matter under 
consideration. For an ultrarelativistic ideal gas this is 

p = cle, (3) 

with Cs = 1/ ^/S being the velocity of sound. For the particular application to heavy-ion 
collisions, we may think of an ideal (massless) pion gas. Then, the energy density is given 
in terms of the temperature T as 

e = g-T\ (4) 
^30 ' ^ ' 

where g = 3 is the degeneracy factor for pions (assuming isospin symmetry). The particle 
density n in the local rest frame is 



n = g^-^T'. (5) 
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In the CM frame, the initial conditions for the one-dimensional expansion of a finite slab 
of length 2 L into vacuum |jlO[ are 




1 



1 



II - ^ (6) 

\z\> L , ^ ' 

— oo < z < — L 

-L<z<L (7) 
L < z < oo . 



Note that all lengths and times scale with L. 

The system evolves as follows, cf. Fig. 1. Initially, the discontinuities at 2; = ±L decay 
into two simple (Riemann) rarefaction waves. The wave fronts move into the slab with 
velocities ±Cs and into the vacuum with the speed of light ±1. At the CM time tm = L/cg 
the two simple waves meet at the center and begin to overlap. At later times t > there 
are three distinct regions of hydrodynamic flow in the forward light cone. We call the one 
at the center of the system, where the two rarefaction waves overlap, the Landau region 



jlOl , |lT| |T2|, because Landau was the first to determine the hydrodynamic solution in this 
region. The other two are called Riemann regions, because here the hydrodynamic solution 
is still given by the original simple Riemann rarefaction waves. 

Since we ultimately want to apply the freeze-out scenarios discussed in this work to 
more complicated, three-dimensional simulations, we study the Landau expansion numer- 



ically. For the numerical solution we employ the explicit SHASTA algorithm The 
applicability of this version of SHASTA to relativistic one-dimensional compression and 
expansion problems in thermodynamically normal matter has been demonstrated in Refs. 



1^ , . We employ a reduced antidiffusion and also a simple full-step method for updating 
source termsQ. 

The initial temperature profile and typical profiles for times t > tm are shown in Fig. 
2. The temperature is calculated from the local energy density e according to @, and 
normalized to the initial temperature Tq. Note that the overlapping rarefaction waves in the 
Landau region produce a local minimum of energy density and temperature at fixed CM 
times t > tm- This is essentially a relativistic effect: faster matter at larger \z\ experiences 
more time dilation and thus cools not as fast in the CM frame as slower matter near the 
center [111 . 



3 The freeze— out scenarios 
3.1 General aspects 

A particle "freezes out" when its interactions with the rest of the system cease. In principle, 
this can happen at any space-time point of the many-particle evolution, depending on the 

■^The half-step treatment of source terms is necessary only to appropriately describe compression shocks 
0- 
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(finite) mean free patli of tlie particular particle. In a one-component system, tlie mean free 
path is estimated as 

A ^ — , (8) 
na 

where a is the (total) cross section and n the (local) comoving particle density. If A is small 
enough such that local thermodynamical equilibrium is established, ideal hydrodynamics is 
valid (and n is, for a pionic system, given by eq. (^). For larger A non-equilibrium effects 
become increasingly important, and if A exceeds the system's dimension D, the latter starts 
to decouple into free-streaming particles. To describe the system's evolution in the latter 
two stages in principle requires kinetic theory. 

Up to now, drastic simplifications of the freeze-out stage were made for practical cal- 



culations [10, 16 1 . The first assumption is that the intermediate region where A ~ Z) has 
negligible extent and can be approximated by a hypersurface in space-time. When a fluid 
element crosses this hypersurface, particles contained in that element freeze out instanta- 
neously. The second assumption is that the freeze-out of particles does not change the 
hydrodynamical evolution. Because of energy and momentum carried away by the leaving 
particles, the hydrodynamical solution is obviously influenced. Also, the decoupling of a 
particle from the flow perturbes (local) thermodynamic and chemical equilibrium in the re- 
spective fluid element. In general, the freeze-out should be treated in a self-consistent way 
taking all these effects into account. This, however, is out of the scope of the present work. 

As mentioned in the Introduction, here we also assume instantaneous freeze-out across 
a hypersurface. We focus on two different scenarios. The first, referred to in the following 
as isothermal freeze-out, corresponds to freeze-out at a constant fluid temperature Ty. The 
second corresponds to freeze-out at a constant CM time tj and is referred to as isochronous 
freeze-out. 

3.2 Isothermal freeze— out 

The isothermal freeze-out is based on assumptions originally formulated by Landau [0 . He 
argued that the freeze-out of an interacting hadronic fluid into free-streaming particles takes 
place at a certain freeze-out temperature Tf which, according to eq. (||) with A ~ D, depends 
only weakly on the system's size, Tf ~ D^^^^. Thus, the isotherm with T = Tf defines the 
freeze-out hypersurface. For pions we have (roughly) a ~ 1 fm^ (at the energy scale relevant 
for pion-pion collisions in a thermal ensemble). We also assume D ~ 2 ~ 20 L ~ 10 fm (cf. 
Fig. 2), where the initial size of the hot central region was estimated to be 2 L ~ 2 to ~ 1 fm 
(for an equilibration time to ~ 0.5 fm after impact), and obtain with eq. (^) the freeze-out 
criterion Tf ~ m^r, where ~ 138 MeV is the pion mass. 

In order to apply this scenario to our underlying expansion model presented in Section 
2, the isotherm of temperature Tf has to be found numerically on the calculational grid. As 
an example. Fig. 1 shows the isotherm Tf = 0.4 Tq. The isotherm is in very good agreement 



with that calculated with the (semi-analytical) method of characteristics, see . 

A hypersurface element is called space-like if its normal vector is time-like, and vice versa. 
The whole isothermal hypersurface can be divided into space- and time-like parts according 
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to Fig. 1. 



3.3 Isochronous freeze— out 

Due to its simplicity, the isochronous freeze-out scenario is widely used in (3+l)-dimensional 
one- and multi-fluid-dynamics The freeze-out condition one commonly assumes is that 
the average temperature of the system at CM time tf drops below a given value Tj, 

ntf) = Tf , (9) 

where T{t) is calculated according to 



m = : , (10) 



where 7 = (1 — t>^)~^/^ and n is the local rest frame density given by eq. 

Note that the hypersurface for isochronous freeze-out is purely space-like. Fig. 1 shows 
the hypersurfaces at CM times tf when (a) the average temperature is T = 0.4 Tq [t ~ 9.2 L, 
lower line), and (b) when the maximum temperature in the system has dropped below Tf 
("late" isochronous freeze-out, t ~ 12.6 L, upper line). Fig. 2 shows the corresponding 
temperature profiles. 

Although the isochronous freeze-out surfaces are purely space-like, we distinguish "time- 
like" (dotted) and "space-like" (full line) parts as shown in Fig. 1. This is motivated by 
the situation at "late" isochronous freeze-out: the energy and momentum associated with 
particles freezing out across in the isothermal freeze-out scenario later has to cross the 
dotted, "time-like" parts of the "late" isochronous freeze-out surface, because it cannot leave 
the forward light cone. On the other hand, the energy and momentum freezing out across 
is the same that freezes out across the "space-like" part of the "late" isochronous freeze-out 
surface. 



4 Calculation of the spectra 
4.1 The Cooper— Frye formula 

In both freeze-out scenarios, the distribution of particles that have decoupled from the fluid 
is determined as follows. The (Lorentz-invariant) momentum-space distribution of particles 
crossing a hypersurface S in Minkowski space is given by |16| 

El^= [ da^k^ fik'^u,) , (11) 



d^k 

where fik^u^) is the local equilibrium distribution function 



/(/c^M^) = — ^ ^ , (12) 
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and do"^ is the normal vector on an (infinitesimal) element of the hypersurface S. do"^ 
is naturally chosen to point outwards with respect to the hotter interior of S, since ([Tl| ) 
is supposed to give the momentum distribution of particles decoupling from the fiuid. In 
our case of one-dimensional geometry, the transverse dimensions enter only as a (constant) 
transverse area factor A. The parametric integration over the hypersurface S in ([TT|) is 
performed by distinguishing space-like parts (with z as integration variable) and time-like 
parts (with t as integration variable), cf. Fig. 1. This yields the total momentum-space 
distribution 

dA^ 
dyk±dk±^ 



where y = Artanh(A;||/£') is the (longitudinal) particle rapidity, m± = (k^ -|- m^)^/^ is the 
transverse mass (which reduces to k± for massless particles), and a = Artanh(i;) is the 
(longitudinal) fiuid rapidity, {dt/dz)^, and {dz/dt)^, are the (local) slope of the space-like 
and the inverse slope of the time-like hypersurface element, respectively. Note that the 
parametric integration in eq. ([T3|) holds for a mathematically positive orientation of S. 



9± 



r ^ smh.y{dt/ dz)Y, — cosh, y 
Jss exp{m^ cosh(a — y)/T} — 1 
sinhy — cosh y{dz/dt)j^ 
St exp{m^ cosh(a; — y)/T} ■ 



+ / dt 



1 



(13) 



4.2 Application to isothermal freeze— out 



We now apply eq. (0) to the isothermal freeze-out scenario. Along the isotherm, T = Tj = 
const, and only a depends on position or time, respectively. The rapidity distribution is 
obtained by integrating ([T3|) over transverse momentum. For massless particles. 



1 dA^ 



K(3) 
27r2 

+ 



To 



dz 



smh.y{dt/dz)j] — coshy 



dt 



cosh^(Q;(z) — y) 
sinhy — cosh y{dz/dt)j^ 
cosh^(a(t) — y) 



(14) 



Analogously, the transverse momentum distribution is obtained integrating over y, 



1 dN 

ATn kidk I 



9 fc± 
47r2 To 



+ 



dt 



dz 



dy 



dy 



sinh y{dt/dz)Y; — coshy 



exp{(/c_L/ro)(ro/T/) cosh{a{z) - y)} - 1 
sinhy — cosh y{dz/dt)Y. 



exp{(fc±/To)(To/T/) cosh(a(t) - y)} 



(15) 



4.3 Application to isochronous freeze— out 

In the isochronous freeze-out scenario the hypersurface S has no time-like part and the 
second term in ( [T3| ) vanishes. Moreover, also the first term in the numerator of the remaining 
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term vanishes due to dt/ dz = Q for this particular hypersurface. The temperature along the 
hypersurface, however, is no longer constant. Thus the rapidity distribution becomes 

^df^^gCi^ rL+tf (T{z)V cosh 7/ 

AT^ dy 27r2 J-L-tf [ Tq ) cosh%a{z) - y) ' ^ ' 

while the transverse momentum distribution is 

ATok^dk^ Att^ To J-L-tf ^J-oo ^ exp{{k^/To){To/T{z)) cosh{a{z) - y)} - 1 ' ^ ' 
Note that the boundaries of the ^-integration correspond to the light cone. 



5 Results 

In this section we discuss the pion rapidity and transverse momentum spectra for isothermal 
and isochronous freeze-out using the expansion model of Section 2. Note that for ideal fluid 
dynamics, entropy is conserved. Since for a gas of massless particles the particle number 
density is proportional to the entropy density, n = 45 C(3)s/2 vr^, also the total particle 
number is conserved. The rapidity distributions presented in this section are numerically 
checked to fulfill this conservation law to better than 2% accuracy (across hypersurfaces 
bounded by the light cone). 

In Fig. 3 we present the rapidity spectrum for isothermal freeze-out, eq. along 
the isotherm Ty = 0.4 Tq that was shown in Fig. 1. The contributions from the different 
parts of the hypersurface as indicated in Fig. 1 are shown separately. The total spectrum 
(full line) shows a double-peak structure, and a small local maximum at y = 0. The time- 
like contributions (dashed) can become negative for particles with a longitudinal velocity 
tanhy which is smaller than the inverse slope of the time-like hypersurface, because then 
sinhy — cosh y{dz/dt)-£ = cosh?/[tanh2/ — {dz/dt)^,] < 0. These contributions correspond 
to particles that are overtaken by the freeze-out condition, and thus in principle re-enter 
the fluid instead of freezing out. Note that it is important to count these particles with 
their correct (negative) sign in order to ensure energy-momentum, entropy, and particle 
number conservation. As can be seen in Fig. 3, the space-like contributions to the spectrum 
(dash-dotted) are always positive, due to \{dt/dz)Y,\ < 1. 

The transverse momentum spectrum, eq. (plSj), is shown in Fig. 4. Here one observes that 
the time-like part of the hypersurface (dashed) contributes less particles than the space- 
like one (dash-dotted), but the slope is the same, reflecting the fact that the temperature 
along the isothermal freeze-out hypersurface is constant by definition. The observed low-A;^ 
enhancement is due to Bose-Einstein statistics. 

For comparison. Fig. 5 shows the rapidity spectrum for isochronous freeze-out, eq. (p!6D, 
with T = 0.4 Tq. The total spectrum (full line) shows a double-peak structure with a local 
minimum ai y = which originates from the time dilation effect discussed in Section 2. 
Different contributions from the "time-like" and "space-like" parts as defined in Section 3.3 
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are shown separately (dashed, dash-dotted hnes). As is obvious from eq. (|TBp, there are no 
negative contributions to the spectrum. 

The corresponding transverse momentum spectrum, eq. ([l^) , is presented in Fig. 6. The 
"space-hke" part of the fluid distribution has on the average a higher temperature and 
thus its contribution to the transverse momentum spectrum a smaller (negative) slope than 
the corresponding spectrum from the "time-like" parts. The total spectrum (full line) is 
nevertheless very similar to the one in the isothermal freeze-out. Fig. 4, because the average 
temperature at freeze-out is the same in both scenarios. 

We now study the dependence of the spectra on the freeze-out temperature. Fig. 7 
shows rapidity spectra for isothermal freeze-out at Tf = 0.8, 0.6, 0.4, 0.2 Tq. For the largest 
T/-values the rapidity spectrum peaks at finite y while for smaller Tf also a peak at central 
rapidity y = develops. Only for Tf -C Tq does this peak dominate the spectrum. 

In the case of isochronous freeze-out the situation is quite different. In Fig. 8, we show 
results for T = 0.8, 0.6, 0.4, 0.2 Tq. For T = 0.8 Tq, one observes a rapidity distribution with 
a Gaussian-like shape, because in this case freeze-out happens so early that the two maxima 
and the minimum at the origin as seen in the temperature distribution (Fig. 2) are too close 
together to be resolved after thermal smearing (especially since in this case the temperature 
is relatively high). Also, for decreasing freeze-out temperatures there are considerable dif- 
ferences to the case of isothermal freeze-out. In particular, there is no maximum at central 
rapidity due to the time dilation effect visible in the temperature distributions. 

We now consider the variation of the transverse momentum spectrum with temperature. 
In Fig. 9 we present the different transverse momentum spectra for isothermal freeze-out and 
for the same freeze-out temperatures as above. As expected, lower freeze-out temperatures 
lead to steeper distributions. The corresponding spectra for isochronous freeze-out are 
shown in Fig. 10. Due to the above mentioned reasons, the results are very similar to those 
of isothermal freeze-out. 

A Gaussian-like shape of the rapidity spectra is often claimed to be a characteristic feature 
of the Landau expansion model |]TU|, |T^. As is demonstrated in Fig. 7, this is not true in 
general and especially not for realistic values of Tf. A Gaussian- like spectrum emerges only if 
one neglects the time-like contributions to the rapidity spectrum, which, however, dominate 
the sum spectra at higher Tf. To illustrate this point we show in Fig. 11 the rapidity 
spectrum originating from the space-like parts of the isothermal freeze-out hypersurface. 
Note that in contrast to Fig. 7, all rapidity spectra have Gaussian-like shapes. The neglected 
time-like contributions are shown separately in Fig. 12. We have normalized the spectra in 
Fig. 11 to their maximum value at ?/ = in order to facilitate a comparison to the results 



of Srivastava et al. [12] who employ this approximation. Note also that their method to 



calculate the spectra violates the conservation laws, since it does not account for energy, 
momentum, entropy, or particle number streaming through the time-like parts of the freeze- 
out hypersurface. 
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6 Summary and Conclusions 



In this paper we have investigated freeze-out in a one-dimensional hydrodynamic model 
for the expansion stage of heavy-ion collisions. We have discussed two different freeze-out 
scenarios used in the literature, namely instantaneous freeze-out across an isothermal and 
an isochronous space-time hypersurface. 

We have found that, provided T|t^ =Tf, the transverse momentum spectra are very simi- 
lar in the two scenarios, but that there are noticeable differences in the rapidity distributions. 
This can be explained as follows. In our model, the fluid evolution is only one-dimensional 
and variations in temperature and fluid velocity occur only along the longitudinal direction. 
It is obvious that they consequently have to vary along different hypersurfaces in space-time. 
These differences are bound to reflect themselves in the (longitudinal) rapidity distribution. 
On the other hand, the transverse momentum spectra average over the longitudinal degrees 
of freedom and thus are only influenced by the (average) temperature. The fact that the 
latter is the same in both scenarios naturally explains the observed similarity. If one, how- 
ever, selects a particular rapidity interval for the calculation of the transverse momentum 
spectra, for instance midrapidity, the main contribution will come from the space-like parts 
of the hypersurface at z=0, where the fluid temperature is different in the two scenarios, 
cf. Fig. 1. The transverse momentum spectra will reflect this difference, cf. Figs. 4 and 6. 
Therefore, one should in general use the (more realistic) isothermal freeze-out scenario in 
three-dimensional calculations. For central collisions with cylindrical symmetry, isothermal 
freeze-out was already employed in the calculations of Ref. . 



We also demonstrated that an improper calculation of the particle rapidity distribution, 
neglecting contributions from time-like parts of the freeze-out hypersurface and thus violat- 
ing the conservation laws, leads to the well-known Gaussian-like shapes commonly associated 
with the Landau expansion model. A consistent calculation, however, produces double- or 
triple-humped distributions. 

Future studies should on the one hand focus on a more realistic description of the de- 
coupling of particles from regions with flnite space-time extent instead across an inflnitely 
narrow space-time hypersurface. On the other hand, the feedback of the decoupling process 
on the evolution of the fluid has to be taken into account to obtain a fully self-consistent 
picture of freeze-out for the hydrodynamical description of heavy-ion collisions. 
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Figure Captions: 



Fig. 1: Illustration of the different freeze-out scenarios. The expansion has been calcu- 
lated with the explicit SHASTA. As an example for the isothermal freeze-out, we show the 
isotherm Tf = 0.4 Tq. It consists of two parts, namely the time-like (dashed), and the 
space-like (dash-dotted). For the isochronous freeze-out, we show the (purely space-like) 
surfaces at the CM times when (a) the average temperature T = 0.4 Tq is reached (lower 
line), and (b) the temperature of the hottest cell has just dropped below Ty ("late" freeze- 
out, upper line). Also for the isochronous scenario we define so-called "time-hke" (dotted) 
and "space-like" (full) contributions, for an explanation see text. The thin dashed lines 
illustrate the boundaries between the Riemann and the Landau region. 

Fig. 2: Temperature profiles for the one- dimensional expansion, normalized to the ini- 
tial temperature Tq. Full line: initial profile. Also shown are profiles at (a) t ~ 9.2 L, when 
the average temperature of the distribution is T = 0.4 Tq (dashed), and (b) at t ~ 12.6 L, 
when the late freeze-out stage is reached (dash-dotted). The profiles have been calculated 
with the explicit SHASTA and standard antidiffusion [Q. 

Fig. 3: Rapidity distribution (full line) for isothermal freeze-out along the isotherm Tf = 
0.4 Tq. The distribution consists of contributions from time-like (dashed) and space-like 
(dash-dotted) parts of the hypersurface. Vq = 2LA is the initial volume of the expanding 
block. 

Fig. 4: Transverse momentum spectrum for isothermal freeze-out at Tf = 0.4 Tq. Time-like 
(dashed) and space-like (dash-dotted) contributions are shown separately. 

Fig. 5: Rapidity distribution (full line) for isochronous freeze-out with T = 0.4 Tq. The 
distribution consists of contributions from the "time-like" (dashed) and "space-like" (dash- 
dotted) parts of the hypersurface. 

Fig. 6: Transverse momentum spectrum for isochronous freeze-out with T = 0.4 Tq. "Time- 
like" (dashed) and "space-like" (dash-dotted) contributions are shown separately. 

Fig. 7: Rapidity distributions for isothermal freeze-out along the isotherms Tf/To = 
0.8, 0.6, 0.4, 0.2. Only the total distributions are shown. 

Fig. 8: Rapidity distributions for isochronous freeze-out at average temperatures T/Tq = 
0.8, 0.6, 0.4, 0.2. Only the total distributions are shown. 

Fig. 9: Transverse momentum spectra for isothermal freeze-out along the isotherms Tf/To = 
0.8, 0.6, 0.4, 0.2. Only the total spectra are shown. 

Fig. 10: Transverse momentum spectra for isochronous freeze-out at average tempera- 
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tures T/To = 0.8, 0.6, 0.4, 0.2. Only the total spectra are shown. 

Fig. 11: Rapidity distributions along the space-hke part of the hypersurfaces for isothermal 
freeze-out with temperatures Tf/T^ = 0.8, 0.6, 0.4, 0.2. The distributions have been nor- 
malized to their maximum value. The time-like parts of the hypersurfaces are neglected. 

Fig. 12: The time-like contributions that have been neglected in Fig. 11. The distribu- 
tions are normalized to the corresponding maxima in Fig. 11. 
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